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ABSTRACT 

The statistical distribution of the magnitude of the vector wind change over 0.25-. 1-, 2-, and 4-h periods 
based on data from October 1995 through March 1996 over central Florida is presented. The wind changes at 
altitudes from 6 to 17 km were measured using the Kennedy Space Center 50-MHz Doppler radar wind profiler. 
Quality controlled profiles were produced every 5 min for 1 12 gates, each representing 150 m in altitude. Gates 
28 through 100 were selected for analysis because of their significance to ascending space launch vehicles. The 
distribution was found to be lognormal. The parameters of the lognormal distribution depend systematically on 
the time interval. This dependence is consistent with the behavior of structure functions in the p tx spectral 
regime. There is a small difference between the 1995 data and the 1996 data, which may represent a weak 
seasonal effect. 


1. Introduction 

This study was motivated by the need to quantify the 
benefit of continuing to operate the Kennedy Space Cen- 
ter (KSC) 50-MHz Doppler radar wind profiler 
(DRWP). In its operational support configuration, the 
instrument measures wind speed and direction at 150-m 
intervals from 2011 to 18 661 m. The data are currently 
used to evaluate wind persistence between the last pre- 
cision wind sounding balloon (“Jimsphere”), about 1-h 
before launch, and launch of a space shuttle. Other 
launch vehicles may soon begin using the DRWP in a 
similar manner. 

The benefit of the DRWP is avoidance of the risk 
associated with undetected wind changes large enough 
to threaten the vehicle or the mission. To assist in as- 
sessing this risk, the statistical distribution of the vector 
magnitudes of 0.25-, 1-, 2-, and 4-h wind changes at 
altitudes between 6 and 17 km was determined over a 
6-month period during the winter season. The 6-17-km 
region is the region in which large launch vehicles are 
most sensitive to wind changes. Winter is the season in 
which the largest short-term temporal changes in winds 
at the target altitudes over Florida are expected. Wind 
changes rather than wind shear were examined because 
space vehicle load programs calculate loads based on 
the wind profiles, and it is the deviation from the ex- 
pected wind profile that poses the hazard regardless of 
whether large shears are involved. 
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Wind profiles were collected every 5 min over a pe- 
riod of 1 17 days from 29 September 1995 through 26 
March 1996 except when the DRWP was down for 
maintenance or repair. The data were extensively quality 
controlled. Vector wind differences were computed at 
each level for the target time intervals (A T) for each 
profile using only data that passed the quality control 
(QC) screening. More than 25 000 profiles were ac- 
cepted for analysis. The first four moments, the prob- 
ability density, and the cumulative probability distri- 
bution of the vector magnitude were computed at each 
level for each A T. 

This paper presents the details of how the study was 
conducted and the results showing the statistical distri- 
butions of the vector wind changes as a function of A T. 
Results include statistics for the distribution at any one 
gate in the target region and also for the maximum 
vector wind change magnitude in the entire region. 

2. The KSC 50-MHz DRWP 

The data used for this study were collected using the 
50-MHz DRWP located near the Space Shuttle Landing 
Facility (SLF). It is a hardware twin to the one at White 
Sands Missile Range described by Nastrom and Eaton 
(1995). The consensus software supplied with the orig- 
inal profiler was replaced by a “first-guess median fil- 
ter” algorithm developed at Marshall Spaceflight Center 
(Wilfong et al. 1993) that also provides capability for 
real-time interactive QC of the data. Details are provided 
in Schumann et al. (1993). 

Profiles are generated every 5 min. Each profile con- 
tains data at each of 112 range gates that represent a 
slice of the atmosphere 150 m thick. Gate 1 is at a 
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nominal altitude of 2011 m while gate 1 12 is at 18 661 
m. The variables measured at each gate include the sig- 
nal power, noise power, Doppler return frequency, and 
spectral width for each of the two oblique beams and 
the vertical beam. These data are converted into esti- 
mates at each gate of the horizontal wind direction and 
speed, the vertical speed, and the signal-to-noise ratio 
in each beam. 

3. The data and its analysis 

a. The data sample 

At each of the 112 range gates, the KSC DRWP pro- 
duces 288 wind speed and direction values daily. Data 
were collected for this study from 29 September 1995 
through 26 March 1996. Excluding days on which the 
DRWP was inoperative, there were 61 days in 1995 and 
56 days in 1996 that yielded 117 days in the sample. 
Some of these days were incomplete due to profiler 
outages, and the QC process removed additional data. 
Since the QC was done gate by gate, the sample size 
varies; but it is typically 25 000 to 28 000 per gate. 
When data are combined for gates 28-100 (6-17 km), 
the resulting sample size approaches two million. 

The time available for this study was limited and a 
full year of data could not be collected or processed. 
The upper air climatology of central Florida is such that 
the winter season brings stronger winds aloft, including 
meandering jet cores with sharp spatial gradients. This 
season was selected for the study as a “worst case” in 
terms of the risk to space launch vehicles. Obviously, 
this is a worst case only in the sense of being worse 
than the summer season. It is unlikely to be the worst 
case over a period of years. 

To estimate the strength of seasonal variation, the 
sample was nearly evenly divided, with the first half 
from 1995 and the second from 1996. Typically, central 
Florida winds aloft and their gradients are somewhat 
stronger in the January through March period than in 
October through December. The mean layer wind speeds 
for gates 28-100 followed the typical pattern. For 29 
September through 31 December 1995 the dataset gave 
a mean speed of 23.19 ms l . For the data from 1 Jan- 
uary through 26 March 1996 the result was 35.46 m 
s '. The respective standard deviations (m s ') were 
1 3.66 and 14.3 1 . The sample size for both sets was near 
a million, and both of the differences are statistically 
significant at the 0.005 level. Comparison of the 1995 
and 1996 results with each other and with the combined 
1995-96 data was used to examine whether a seasonal 
effect was significant during these months. 

b. The quality control (QC) process 

Each range gate in each profile constitutes one record 
in the daily data files. Each record includes a QC flag 
to track individual quality control tests. The active flag 
bits are listed below in Table 1. 


Tabi.k 1. Quality control flag thresholds. 


What causes flag 

Threshold 

value 

Excessive shear internally computed by 
DRWP 

0.1 s 1 

Excessive vertical speed 

3 m s 1 

Excessive spectral width 

3 m s 1 

Exceeds first-guess propagation limit 

6 

Missing profile — data place holders (999) 

N/A 

Fails small median test 

Sec text 

Excessive adjacent gate directional shear 

30° per gate 

Excessive adjacent gate speed shear 

7ms 1 per gate 

Insufficient signaLto-noise ratio (SNR) 

— 15 dB 

Fails manual QC 

See text 


The vertical speed threshold is used to identify pre- 
cipitation. Heavy rain contaminates the radar signal and 
produces bad winds. Typically, heavy rain produces a 
large vertical wind indication based on the radar track- 
ing the falling raindrops. 

The shear value computed internally by the radar soft- 
ware provides a coarse warning that the measurements 
are invalid and would usually not be triggered except 
by an aircraft in the beam or a major failure in the 
system. The more stringent speed and direction shear 
tests are designed to eliminate meteorologically un- 
realistic winds produced by sidelobe or interference sig- 
nals and to flag those signals. 

The small median test is the one described in detail 
by Carr et al. (1995). It examines not only continuity 
of the winds in the vertical but temporal continuity as 
well. The test requires threshold values be designated 
for three different heights. This work used 5.7 m s 1 at 
2 km, 10.2 m s 1 at 9 km, and 8.4 m s 1 at 16 km based 
on about 20 years of local windsonde data. These criteria 
are more restrictive than those used by Carr et al. ( 1 995). 

On many of the days during the experimental period, 
personnel of the Applied Meteorology Unit (AMU) ex- 
amined the spectra in real time and made notes of any 
obvious sidelobe or interference signals that could ad- 
versely affect the data. When the Titan program con- 
ducted real or simulated launch operations, the data were 
not only monitored but were manually quality controlled 
before transmission to the data system. In these cases, 
the real-time manual QC operator made extensive notes. 
Both the AMU notes and the Titan notes were used to 
assist the after-the-fact manual QC process described 
next. 

After the automated QC was run, the file was dis- 
played graphically. Wind speed and wind direction were 
separately but simultaneously presented as time-height 
diagrams in which each range gate at each 5-min in- 
terval was plotted as a colored point. The resulting color 
contour map permitted visual detection of wind speed 
or wind direction changes and the identification of nat- 
ural features like jet streaks and artifacts like sidelobes. 

AMU and Titan QC logs were used for guidance. 
During visual inspection of the data, the time and al- 
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Fig. 1. The cumulative probability (vertical axis) of the magnitude 
of the 1-h wind change at any level in the region from gate 28 to 
gate 100 inclusive exceeding the specified (horizontal axis) value. 
Probabilities are expressed as normalized standard deviations z (e.g., 
a probability of 0.5 equates to ~ - 0. 0.8413 equates to c = 1 )• The 
points are the measured cumulative distribution. The line is the log- 
normal distribution derived from the moments as described in the 
text. 

titude of any apparent sidelobe or interference signatures 
were noted. Data that appeared to be unreliable were 
manually flagged by setting a QC flag bit. 

Nearly all of the manual QC consisted of flagging the 
interior (in time-height space) of sidelobe and interfer- 
ence signatures whose boundaries were flagged by the 
automated process. Of the 117 days in the sample, 44 
required some manual flagging. Less than 1% of those 
data were manually flagged. 

Once a data file had been accepted, the delta files 
were generated and examined visually. Delta files are 
identical to the data files except that they contain the 
wind changes rather than the winds. Sometimes a side- 
lobe or interference signal stood out better in the wind 
changes than in the winds themselves. When this oc- 
curred, the process described above was repeated. When 
the data file had been reaccepted, the delta files w F ere 
rerun and reexamined. The data were used for analysis 
only when the basic data and the delta files had been 
examined and accepted. During analysis, any record for 
which the QC flag is nonzero was ignored. 

c. The analysis scheme 

The initial experimental design called for computing 
the first four statistical moments (mean, variance, skew- 
ness, and kurtosis) and the probability density and dis- 
tribution of the magnitude of the vector wind change at 
each level in the region of interest. It became quickly 
obvious that this was impractical. Greater sample sizes 
with all of their advantages for statistical analysis could 
be obtained by combining data from multiple levels while 
simultaneously reducing the information output to a man- 
ageable size. Limited resources made this strategy nec- 


essary despite the loss of information about the variability 
of the results with height. The height levels from gate 
28 to gate 100 were selected for their relevance to as- 
cending launch vehicles. In addition, the maximum value 
for the magnitude of the vector wind change anywhere 
in the target region was computed for each profile. The 
same statistics were computed for this dataset. 

A combination of gates 32-34 (a shuttle’s transonic 
region) was created to permit "quick-look” examination 
of the extent to which vertical inhomogeneity might be 
a concern. The mean wind speed in this layer averaged 
30% lower than in the layer from gate 28 to gate 100. 
The wind difference data from the two regions differed 
by an amount that was statistically significant due to the 
enormous sample sizes involved, but the differences ap- 
peared insignificant to the author from the point of view 
of risk assessment. After determining that the distribution 
of these data had roughly the same shape and parameter 
values as those in the thicker layer, detailed analysis of 
this dataset was deferred due to limited resources. Lim- 
ited resources also prevented any attempt to stratify the 
data temporally to look for diurnal effects. 

Examination of the statistical results from the first two 
days of the September 1995 data suggested that the data 
were distributed lognormally. The experimental design 
was enhanced to include generation of six estimates of 
the parameters of a lognormal distribution from the four 
moments by combining them in pairs. Each pair of mo- 
ments generates an estimate for the two lognormal pa- 
rameters as shown below. The mutual consistency of 
these estimates served as an indicator that the distribu- 
tions really were lognormal. This was confirmed by using 
the mean of the six estimates of each parameter to gen- 
erate a straight line on a log probability plot on which 
the actual distribution was overlaid for comparison. 

The formula for the nth central moment of a lognormal 
distribution is given by Aitchison and Brown (1966) as 

M n — exp(nyit + ncr 1 ! 2), (1) 

where jjl is the mean of the logarithms and cr is the 
standard deviation of the logarithms of the lognormal 
variable. For any two values of n for which M n are 
known, the result is a pair of simultaneous linear al- 
gebraic equations for fx and <r 2 , which may be solved 
readily. 

4. The results 

This section presents the analysis showing that the 
magnitude of the vector wind change is lognormally 
distributed. The variation of the parameters of the log- 
normal distribution with the length of the time interval 
over which the change takes place is also presented. A 
brief discussion of possible seasonal variation in the 
distribution is included. Finally, the implications of 
these results for spaceflight operations and for the design 
of profiler networks are discussed. 
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Fig. 2. The cumulative probability of the largest single value of 
the l-h wind change in the region from gate 28 to gate 100 inclusive 
exceeding the specified value. The format is the same as in Fig. 1. 


a. The magnitude of the vector wind change is 
lognormally distributed 

Figure 1 shows the cumulative distribution of the 
combined gates 28-100 data for the entire sample of 
l-h wind changes. The individually plotted points are 
the actual measured distribution. The straight line rep- 
resents a pure lognormal distribution with parameters 
equal to the average of those derived from the six pair- 
wise combinations of the first four moments of the mea- 
sured distribution. The horizontal axis is the magnitude 
of the vector wind change over a 1 -h period. This sample 


contained 1 882 502 data points. To the extent that the 
wind differences may be related to wind variances, these 
results agree with those of Nastrom and Eaton (1995), 
who found the variances over a l-h period at similar 
altitudes were lognormally distributed. 

Figure 2 shows the same information as Fig. 1 except 
that the data are the maximum values of the magnitude 
of the vector wind change in gates 28-100 for each 
profile rather than the individual gate values. The sample 
size is 25 983. This figure shows that the maximum 
values in the entire region are also lognormally distrib- 
uted. 

The data from 1995 and 1996 separately and the in- 
dividual monthly data behave in the same manner, as 
do the data for time periods of 0.25, 2, and 4 h. Figures 
1 and 2 above are typical of any of these data subsets. 

b. The parameters of the lognormal distribution 

depend on AT 

The two parameters of the lognormal distribution, p 
and a ; were computed for velocity magnitudes in meters 
per second. Figure 3 shows the variation of p with the 
time over which the changes are measured (AT) for the 
gates 28-100 combined data. Figure 4 shows the vari- 
ation of p for the layer maxima. Despite the significant 
scatter, a clear, roughly linear relation between the value 
of p and the logarithm of AT is apparent in both cases. 
Linear regressions of p versus ln(AT) yield r > 0.9 
(significant to <0.005) in all cases. 

Figures 5 and 6 present the same information for a. 



Fig. 3. The relationship between the mean ( fi ) of the natural logarithm of the magnitude of 
the vector wind change (ms 1 ) change at any height level in the region from gate 28 to gate 
100 inclusive and the time interval (A T) over which the change takes place. The lines represent 
combined data from September-December 1995, January-March 1996, and the entire dataset. 
The points represent individual months to show the scatter. 
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Flg. 4. The relationship between the mean </x) of the natural logarithm of the magnitude of the 
largest single vector wind change (m s 1 ) change in the region trom gate 28 to gate 100 inclusive 
and the time interval (AD over which the change takes place. The format is the same as in Fig. 3. 


Although the scatter is much larger, a similar relation 
with opposite sign is suggested. Linear regressions of 
c r versus ln(A7) yield r 2 > 0.4 (also significant to 
<0.005) in all cases. 

The relationships between the lognormal parameters 
and AT can be used to interpolate these results to any 


AT between 0.25 and 4 h. They can also be used with 
appropriate caution to extrapolate these results to larger 
and smaller time intervals. Extrapolation to larger time 
intervals is less justified because semidiurnal, diurnal, 
and synoptic features should begin to dominate sto- 
chastic ones at these scales. 



Flg. 5. The relationship between the standard deviation (<r) of the natural logarithm of the 
magnitude of the vector wind change (ms ') change at any height level in the region from gate 
28 to gate 100 inclusive and the time interval (AT) over which the change takes place. The 
format is the same as in Fig. 3. 
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c. The implications of the distribution 

Examination of the previous figures suggests that part 
of the scatter may be due to secular changes in the 
environment during the experimental period. The p and 
a values tend to change in opposite directions with A T. 
The p values for 1996 were slightly higher than for 
1995, but the cr values were lower, thus maintaining the 
same negative correlation. Although the scatter is large, 
a t test shows the differences in p are statistically sig- 
nificant at the 0.001 level, and a chi-squared test shows 
that the differences in crare significant at the 0.005 level. 
In any case, the secular variation is not large compared 
with the scatter and it remains to be seen whether it is 
meteorologically significant. These effects, whatever 
their cause, also naturally appear in the probability of 
exceedance curves presented in the next section. 

The lognormal distribution of the magnitude of the 
vector wind changes facilitates comparison of these re- 
sults with atmospheric structure functions. Structure 
functions for wind are usually defined for each velocity 
component in terms of spatial separation by or in terms 
of temporal separation by 

DJr) = ([u,(x) - w ,( a + r)] 2 ) (2) 

or in terms of temporal separation by 

D n (r) = ([14,(1) ~ u,(t + r)] 2 ), (3) 

where the brackets denote the average and the subscript 
i denotes the #th component (Stull 1989). 

If the logarithm of the wind differences is normally 
distributed, then the logarithm of their square is also 


normally distributed with twice the mean and twice the 
standard deviation of the distribution of the log of the 
differences. If the parameters of the distribution are 
known as a function of temporal separation, then the 
mean value of the squared differences can be calculated 
as a function of the temporal separation. Since the whole 
vector magnitude was used in this study rather than any 
single component, the result will be a “magnitude struc- 
ture function’' equal to the sum of the three component 
structure functions. 

The small number of independent estimates of the 
behavior of the distribution parameters with temporal 
separation and the scatter in the data suggest that pre- 
cision calculations are unwarranted. To compare these 
results with those published in the literature, the vari- 
ation of p and cr with AT was approximated by a straight 
line fit through the 0.25- and 4-h values for the com- 
bined dataset. This resulted in linear equations for p 
and a as follows: 

2p = a M + b M InA / (4) 

2cr = a (T + b ir lnAr, (5) 

where the factor of 2 converts the fit for differences to 
one for squares of differences. Time differences in sec- 
onds were used for dimensional consistency. 

The resulting mean value of the squared differences, 
D(A/), may be reduced to the form 

D(At) = C,Ar v A/ a h,A/l , (6) 

where C t = exp(r/ M 4- 0.5r/ 2 ), N = b^ - a t Jb in and k = 
0.5fr 2 . When the data from Figs. 3 and 5 are used to 
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Fig. 7. The probability (vertical axis) of the magnitude of the 1-h wind change at any height 
level in the region from gate 28 to gate 100 inclusive exceeding the specified (horizontal axis) 
value for various values of AT. The lines are results from the entire dataset. The points separately 
present the 1995 and 1996 data to show the variation from late tall to early spring. 


generate the constants in the linear equations given 
above, the result is C, = 0.01, N — 0.8, and k = 0.001 
to one significant figure. 

Over the range from 0.25- to 4-h, the term containing 
k varies from near unity to 1.17. If we treat this term 
as approximately constant, the result resembles that 
widely reported in the literature for structure functions 
in the inertial subrange if the assumption is made that 
spatial and temporal variations may be treated similarly 
using a Taylor's hypothesis. The spatial structure func- 
tion under these conditions is given by 

Dir) = Ce m r v \ (7) 

where C is a constant of order unity (Ottersten 1969). 
The value of N obtained here is within 20% of the 
inertial subrange value. Moreover, the value of C, ob- 
tained here implies that e is of the order of 0.001 m : 
s \ which compares well with values obtained in the 
free atmosphere from aircraft observations (Merceret 
1976; Vinnichenko and Dutton 1969), although up to 
two orders of magnitude smaller values have been re- 
ported in the upper troposphere and lower stratosphere 
(Vinnichenko 1970). Significantly, the data obtained 
here were taken during the winter where both the wind 
speed and wind shear are greatest, suggesting the largest 
degree of production and dissipation of mechanically 
generated turbulence. 

No spectral analysis has yet been undertaken with 
these data, but Vinnichenko (1970) concluded that in 
the mesoscale (periods from 1 day to about 10 min), 


the —5/3 power law is followed fairly well. The data 
here fall within those scales, and the slope of the struc- 
ture function is consistent with that spectral behavior. 

d. Significance for spaceflight operations and the 

design of upper -air networks 

Using the lognormal model with parameter values 
based on the data collected in this study, we have gen- 
erated the required curves for 0.25-, 1-, and 4-h inter- 
vals. The data for the 2-h interval were also generated, 
but are not plotted in the figures to enhance legibility. 
Figure 7 presents the curves for probability of exceed- 
ance at any preselected gate from 28 to 100. The vertical 
axis is the probability that the magnitude of the vector 
wind change over the specified time period will exceed 
the corresponding value on the horizontal axis. Figure 
8 presents the probability that the threshold will be ex- 
ceeded anywhere in the entire layer. 

For spaceflight operations, the benefit of rapid access 
to wind profile data is the avoidance of unanticipated 
wind changes, which pose the threat of destruction of 
a launch vehicle (Mecham 1995). For assessing the haz- 
ards to space launch operations, engineers must first 
determine what wind change between the last balloon 
and launch poses an operational threat to the launch 
vehicle. The curves presented here can then be used to 
quantify the likelihood of that threat. If the risk is un- 
acceptably high, then availability of day of launch 
DRWP capability should be a mission requirement. A 
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Fig. 8. The probability (vertical axis) of the magnitude of the largest l-h wind change in the 
region from gate 28 to gate 100 inclusive exceeding the specified (horizontal axis) value for various 
values of IT. The format is the same as in Fig. 7. 


larger dataset and additional analysis could be used to 
stratify the risk by altitude and season. 

For design of upper-air networks, the required tem- 
poral sampling rate depends on the need to assume that 
significant features are observed and the data are rep- 
resentative (Douglas and Stensrud 1996). Sampling 
rates for current profiler networks are typically one per 
hour based on consensus averages of 10 or more profiles 
(NOAA 1994; Spencer et al. 1996). These results allow 
quantitative assessment of the wind changes that might 
be missed with sampling or reporting intervals greater 
than a specified time. This information can be used to 
optimize the trade-off between cost and complexity of 
data collection and the accuracy and completeness of 
the data. 

5. Summary and conclusions 

An examination of six months of winter season 
50-MHz DRWP data over central Florida has shown 
that the magnitude of vector wind changes for periods 
from 0.25 to 4 h is lognormally distributed. The log- 
normal parameters vary systematically and simply with 
the time interval. 

Probability of exceedance curves for wind changes 
at any single 150-m gate in the region between 6 and 
17 km and also for the maxima within the entire region 
were also generated. These curves can be used to assess 
the risk to space launch programs of short period wind 
changes and to assist in determining the required tem- 
poral resolution for proposed profiler networks. They 


can also be used to assist in the determination of w hether 
the cost of operating a profiler is warranted for a par- 
ticular application given a specified acceptable risk 
level. 
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